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The Polyakov loop and the chiral condensate are used as order parameters to explore analytically the possible 
phase structure of finite temperature QCD. Nambu-Jona-Lasinio models in a background temporal gauge field are 
combined with a Polyakov loop potential in a form suitable for both the lattice and the continuum. Three possible 
behaviors are found: a first-order transition, a second-order transition, and a region with both transitions. 



1. INTRODUCTION 

Our recent workjlj on an effective action for 
dynamical quarks has caused us to examine the 
simple picture in which dcconfincment and chiral 
symmetry are independent. On the basis of this 
work, it appears necessary to consider models of 
the deconfining and chiral symmetry transitions 
in which the two order parameters, the chiral con- 
densate and the Polyakov loop, are coupled. 

In this paper we study a two flavor Nambu- 
Jona-Lasinio model p in a uniform background 
temporal gauge field with a Euclidean action of 
the form 

C = + m)# + G [(vj/*) 2 + (vj/i 75 f*) 2 ] (1) 
where 



dn+igA^x) with A^(x) = S 0fi A . (2) 



The spontaneous breakdown of Z$ symmetry in- 
duced by the gauge field is modeled by a simple 
potential with cubic and quartic terms. Com- 
bined with the quark interaction effects, this 
model explicitly constructs a free energy depen- 
dent on the two order parameters, in the spirit of 
Landau-Ginsberg arguments. Because this model 
generates an effective potential consistent with 
two- flavor QCD, we will show that the class of 
critical behaviors possible in finite temperature 
QCD is a priori larger than has been previously 
realized. 
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2. EFFECTIVE POTENTIAL FOR FI- 
NITE TEMPERATURE NJL MODEL 

We let S denote the propagator for a free 
fcrmion with current mass m and S denote the 
propagator for a fermion with constituent mass 
M. Following Cornwall, Jackiw, and Tomboulis 
||, a one-loop effective action r(5) for the theory 
given in Eq. (|l|) at finite temperature is given by 

r(5) = Tr [In (S^S)] - Tr (S^S - l) 



G 



dt / d^xlTrSix,!)}' 



(3) 



where G = [1 + l/(2N f N c )]G. The effective po- 
tential can be conveniently written in terms of the 
constituent mass M: 
(M - m) 2 



V(M,T) = 



2 J1l v 

kn 



AG 
d 3 k 



-tr r In 



(ko+gAof+Lol (4) 



where lot — + M 2 . Evaluating the mode 
sum over ko, we find that up to an irrelevant con- 
stant 



4G 



(2tt)3^ 



2N f f d 3 k 




2N S 



(2tt) 3 
d 3 k 



ti' c In (l 



-0v, 



V) 



J (2tt)! 



tr c In (1 + e^^) 



(5) 
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where V = e l ^ aA ° is the Polyakov loop associated 
with the background A . 

We regulate V(M) by introducing a non- 
covariant cut-off, 9 (A 2 — k 2 ), and make the ap- 
proximation 



tr c (7>") « N c 



tT c V 



(G) 



Noting that the quark determinant forces the 
trace of the Polyakov loop to be real, Eq. (||) now 
becomes 



V(M, T) 



2N f N c 

1T 2 (3 



(M - m) 2 NfN c 

AG 7T 2 



dfc fc 2 In 



i-/3u> k 



dk k u>k 



tr c V 



(7) 



We determine A, M, and G by simultaneously 
fixing the chiral condensate and the pion form 
factor (|, with < ^ >= N f (246.7 MeVf and 
/ w = 93MeV. For the remainder of this paper we 
will use the numerical solution 



M = Hh.lMeV 



GA = 2.196, 



A = 630.9MeV 



(8) 



which assumes a current quark mass m u — rrid = 
AMeV. 

3. POLYAKOV LOOP POTENTIAL 

We denote tr c V by (p. Near the deconfinement 
temperature To the potential of Polyakov loops in 
pure QCD can be approximated by a polynomial 
potential of the form flfl 



u{<t>, t) = aU 



T-T, 



D 



TnL- 



>D 



(9) 



The potential is parametrized such that for a pure 
gauge theory the phase transition occurs at To 
with a latent heat L and the Polyakov loop jump- 
ing from to 4>d- 

This choice for U is phenomenological. Al- 
though perturbative QCD does yield a quartic 
polynomial potential for the Aq field, the poten- 
tial so obtained does not display critical behavior, 



and it is only valid at high temperatures . How- 
ever, it does suggest that the parameter A should 
be taken to be on the order of T D . We have chosen 



A = T% and L 



2T D for the numerical calcula- 



tions below, neglecting any possible temperature 
corrections or other dependencies in A, Tjj, 4>d 
and L. The direct association of To, 4>d and L 
with measureable quantities holds only for very 
heavy quarks. For light quarks, these parameters 
must be determined by fitting to the observed be- 
havior. 

We can now couple the chiral symmetry and 
deconfinement phase transitions. Let 

V(M, 0, T) = V(M, & T) + U(<f>, T), (10) 

where V(M, <j>, T) is given in Eq. (Q) with tr c V = 
<f). To determine the critical behavior of the cou- 
pled effective potential, the absolute minimum of 
the potential V(M, (j>, T) as a function of M and 
<f> must be found as T is varied. A satisfactory de- 
termination of the entire phase diagram requires 
numerical investigation. 

4. RESULTS 

The two flavor Nambu-Jona-Lasinio model has 
a second-order chiral phase transition for mass- 
less quarks. For the parameter set used here, this 
transition occurs at T x = 194.6 MeV, as deter- 
mined numerically from V(M,<f> = 1,T). The 
order of the transition is consistent with Monte 
Carlo simulations of two- flavor QCD, and the 
critical temperature is plausible. However, the 
finite temperature quark determinant is not con- 
sistent with the known behavior of QCD unless 
the effects of a non-trivial Polyakov loop are in- 
cluded. As Eq. U makes clear, the effects of fi- 
nite temperature in the quark determinant are 
suppressed by the small expected value of the 
Polyakov loop at low temperatures. Without the 
Polyakov loop effects, the conventional Nambu- 
Jona-Lasinio model displays the T 4 behavior of a 
free quark gas at arbitrarily low temperatures. 

Numerical investigation shows that the critical 
behavior of V(M, 0, T) is most sensitive to To and 
(j)o ■ Figure |l| illustrates the three distinct types 
of critical behavior observed as To and <f>o are 
varied with the other parameters held fixed and 
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merit, and a later second-order chiral symmetry 
restoring transition. In Region III the Polyakov 
loop experiences a first-order jump at some T c , 
and simultaneously the effective quark mass drops 
suddenly to MeV, indicating a single first-order 
transition which restores chiral symmetry. 

Figure || shows the location of the three regions 
in the (To, </>d) plane. Figure || plots the behav- 
ior of the Polyakov loop, <f), and the constituent 
quark mass, M, as a function of temperature at 
the point (T D = 220MeV, <j> D = 1.2) in Region II. 
There is a clear separation of the first and second- 
order transitions. 
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Figure 1. Schematic behavior of <j> and M in the 
three regions of the (Td,0£>) plane as tempera- 
ture increases. 



the current mass m set to 0. In Region I of the 
(To, plane the effective quark mass changes 
continuously from its T = value to a mass of 
MeV as the temperature increases, signaling a 
second-order chiral phase transition. Simultane- 
ously, the expectation value of the trace of the 
Polyakov loop increases smoothly with tempera- 
ture, exhibiting a large, but continuous, increase 
in the vicinity of the chiral transition. This region 
exhibits behavior most similar to that observed in 
simulations of two- flavor QCD. In Region II the 
Polyakov loop experiences a first-order jump at 
some T c . At the same T c the constituent quark 
mass undergoes a sudden drop, but chiral sym- 
metry is not restored. As the temperature in- 
creases beyond T c , the constituent mass moves 
continuously to zero. Thus, there are two dis- 
tinct phase transitions in Region II, a first-order 
transition driven by the dynamics of deconfine- 
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Figure 2. Phase diagram in the (4>d,Td) plane. 



5. LATTICE VERSION OF THE EFFEC- 
TIVE POTENTIAL 

A lattice version of the effective potential has 
the advantage that many parameters can be set 
from QCD simulations. There are two related 
issues in extending our work to the lattice: the 
choice of a lattice fermion formalism and setting 
other parameters in the effective action. 

A lattice version of the NJL model can 
be straightforwardly implemented using naive 
fermions || . Species doubling can be handled by 
a formal replacement of N c with N c /16. Using 
the single quark current mass, chiral condensate, 
and constituent mass of the last section as inputs, 
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Figure 3. Plot of <j> and M as a function of tem- 
perature at the point (T D = 220MeV, (j) D = 1.2) 
in Region II. 



we find 

a- 1 = 389.0MeV Ga~ 2 = 0.8343. 



(11) 



Hence, f n = 94.09 MeV. Also note that M T=0 a 
is of order 1. The temperature is given by 
T = l/(N t a) where N t is the temporal extent 
of the lattice. This places a problematic upper 
limit on a discrete range of available tempera- 
tures. Given our choice of parameters, T max is 
only 194.5 MeV. 

The temperature may be changed continuously 
by an asymmetric rescaling of the lattice in which 
the lattice spacing at in the temporal direction 
varies independently of the spatial lattice spac- 
ing a s ||. Defining a s = a and £ = cit/a s , we 
now have in physical units T — l/(N t a^). The 
effective potential is 



V(M,<t>,0 




N f N c I 
4£ J. 


* d 3 k 
-« (2^) 3 


NfN c 1 

J. 


ff d 3 k 


- 1 n 


tr c {4>) -i 



In 



l + £ 2 A 2 (k)+£A(k) 



l + e 2 A 2 (k) 



4G 



(12) 



where A(k) = J M 2 + J2j sm (Pj)- However, us- 
ing the same inputs as the symmetric lattice case, 



as well as /„• = 93 MeV, we find the unique an- 
swer 



477.2MeV 



Ga 



1.255 



(13) 



with £ = 2.002. 

An alternative to asymmetric lattices might 
be variant or improved actions for which a" 1 is 
larger. 

6. CONCLUSION 

Polyakov loop effects have a strong impact on 
the possible critical behavior of the NJL model. 
The universality argument [jl0| which predicts a 
second-order chiral transition for two flavors and 
a first-order transition for three or more flavors 
may fail when this additional order parameter is 
included. 
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